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As a first step towards the understanding of the thermody- 
namical properties of quasiperiodic structures, we have per- 
formed both analytical and numerical calculations of the spe- 
cific heats associated with successive hierarchical approxima- 
tions to multiscale fractal energy spectra. We show that, in a 
certain range of temperatures, the specific heat displays log- 
periodic oscillations as a function of the temperature. We ex- 
hibit scaling arguments that allow for relating the mean value 
as well as the amplitude and the period of the oscillations to 
the characteristic scales of the spectrum. 
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I. INTRODUCTION 

The discovery of quasicrystals in 1984 [Q aroused a 
great interest in quasiperiodic structures, as is confirmed 
by the great number of theoretical B and experimen- 
tal j3| works that followed (see also ||]). In particular, 
the behavior of a variety of particles and quasi-particles 
(electrons ||, phonons and others Q) in quasiperi- 
odic structures has been and is currently being studied. 
A fractal energy spectrum is a common feature to such 
structures (e.g., ||). As in general these spectra tend 
to be quite complex, simple models have been studied 
to enlighten the thermodynamical specificities that such 
systems may display. Within this vein, we studied in |^] 
one of the simplest fractal spectra (the triadic Cantor 
set); there it was shown that the specific heat of such 
system exhibits a very particular behaviour: it oscillates 
log-periodically around a mean value which equals the 
fractal dimension of the spectrum. 

In this paper we extend the analysis of Q to the case 
of multi-scale spectra and present a theoretical connec- 
tion between the structure of the spectrum and the cor- 
responding thermodynamical properties. In particular, 
we make explicit the relationship between discrete scale- 
invariance and log-periodic oscillations. We start our 
discussion with the two scale (r±, r2) fractal set, in its dis- 
crete version (Section II). Then, in Section III, we show 
that the results of II can be extended to the continuous 
two-scale case, as well as to the n-scale problem. Section 
IV contains the conclusions. 



II. TWO SCALE DISCRETE MODEL 

Let us begin by considering a spectrum lying on the 
{tx,t%) two scale fractal set |T^] indicated in Fig. la (the 
restrictions < r\ , ri , r\ + T2 < 1 apply) . The starting 
point (n = 0) for the construction of this spectrum is an 
arbitrary discrete set of levels lying in an energy interval 
we take to be [0,1]. This is the pattern that will be 
repeated ad infinitum in a self-similar way. 
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FIG. 1. Energy spectra. The first four steps in the con- 
struction of (ri,r2) fractal sets. We show a discrete case (a) 
in which the starting point (n = 0) is a set of two levels at 
e = 0, 1. These levels are then compressed by a factor r± (r2) 
and put on bottom (top) of the interval [0, 1] (n = 1), and so 
on for increasing n. The construction of a continuous exam- 
ple (b) starts from a band of uniform density in [0, 1]. The 
iterative rule is the same as in (a), i.e. n — 1 corresponds 
to a spectrum whose first and second bands are the intervals 
[0, n] and [1 — ra, 1] respectively, etc. We take the level den- 
sity inside each band to be a constant, and the same for all 
bands in a given hierarchy. In both cases (a) and (b) , a fractal 
emerges at the n — > oo limit. 

The next step consists in compressing the set n = 
by factors n and r-i and putting the two resulting pieces 
on bottom and top of the interval [0, 1], respectively (see 
Fig. la). Recursive application of this rule eventually 
leads to a set of fractal dimension df given by r 1 / +r 2 f — 
1 (hence, if r\ = T2 = r, df = — In 2/ In r). This rule 
can be explicitly written as a recurrence equation for the 
energy levels e,- at the n-th stage, 
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This analytical rule for the construction of the spectrum 
is the key to obtaining scaling relations for the thermo- 
dynamical quantities. The starting point is the partition 
function for a given hierarchy n: 
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where (3 is the inverse temperature (we are considering a 
unit Boltzmann constant, i.e., kg = 1). A normalization 
prefactor has been included so that is well defined 
in the limit n — > oo; in any case, it will not affect the 
thermodynamics of the system. Now, a recurrence for- 
mula for the partition function is readily obtained as a 
direct consequence of the self-similarity of the energy set 

Z (n+i)( /3 ) = \zW(pn) +e- /3(1 ~ r2) Z(™)(/3r 2 )] /2 . (3) 
Introducing Z(/3) = lim„^oo Z^ n \(3), we have 

Z((3) = [Z((3n) + e-^-^Ztfrz)] /2 . (4) 

From now on we will restrict our discussion to the low 
temperature regime, as this is the most interesting one. 
In fact, as the temperature is lowered, the smaller scales 
of the fractal are progressively revealed, and anomalous 
effects are expected. Moreover, in the case T « 1 the 
analysis gets simplified and some general conclusions can 
be obtained. In this regime we can safely neglect the ex- 
ponentially small term in (|j) and derive scaling relations 
for the partition function, the dimensionless free energy 
Q = F/T = —\nZ, the total energy E, the entropy 5, 
and the specific heat C: 



Z{T) = Z(T/n)/2 
Q(T) = Q(T/n) + In 2 
E(T) = ri E(T/n) 
S(T) = S(T/n) - In 2 
C(T) = C(T/ n ) . 



(5) 
(6) 
(7) 
(8) 
(9) 



Independently of the n — energy pattern, the relevant 
scale factor is r\ (conversely, r 2 governs the scaling laws 
for negative temperatures). 

The most interesting of the equalities above is the last 
one, which expresses the fact that the specific heat is a 
log-periodic function of the temperature, that is C(T) — 
/(27rlnT/ lnri), where / is a 2-7r-periodic function. In 
other words, if one sets T — rf, C(x) results in a periodic 
function of x (of period one). Consistently, its mean value 
can be calculated as 



/•So+l 

(C(T))= C(rf)dx-. 

Jxq 

= S(t) - Sjrrr) 
In r\ 
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In Ref. 1 it was shown that for a one-scale Cantor spec- 
trum (i.e. n = r 2 = 1/3), the average value (C(T)) 
coincides with the fractal dimension df. The equality 
above shows that this is not the case for a two scale frac- 
tal. We will come back to this point later on to argue that 
the "dimension" d = — In 2/ lnr i can be given a simple 
meaning. We remark that the result (|l0|) holds as long 
as T -C 1 (we recall that, as the spectrum is bounded, for 
high temperatures the specific heat must decay as T~ 2 ). 

The scaling reasoning has given us information con- 
cerning the mean values. In order to discuss the oscil- 
lations around the mean value, we resort to a numerical 



analysis. Starting from (f 
proximations to C(T), 

d_ 

df 



we have computed finite ap- 
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and studied its dependences on the hierarchical depth n, 
and the parameters n and r%. In Fig. 2 we show some 
plots of the specific heat vs. temperature for different 
values of (ri,r 2 ) and fixed hierarchical depth n = 8. 




FIG. 2. Specific heat (in units of ks) vs. temperature (in 
units of the width of the spectrum) for the (ri,ra) discrete 
fractal set of Fig. 1 (n = 8). Two levels at e = 0, 1 were 
taken as the n = pattern. The curves are parametrized 
by the scale factors (ri, r%). The horizontal lines indicate the 
average value (C) = d — — In 2/ Inn. The dotted lines corre- 
spond to our prediction C w d+a" cos(o; lnT) + 6" sin(oi InT), 
where ui = — 27r/lnri. The parameters a", b" are related to 
basic properties of the smoothed spectrum (see text). For 
high temperatures (InT > 0) the specific heat decays as T -2 , 
for arbitrary n. The low-temperature breakdown of the oscil- 
latory behaviour is pushed towards the left when n increases. 



It is 



apparent that there is range of temperatures in 
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which C(T) behaves in the way the above scaling ar- 
guments predict. This is a range of intermediate tem- 
peratures, T m i n <C T <C 1 where T mm ~ r" is asso- 
ciated to the smallest scale of the "truncated fractal"; 
of course, T mnl — > in the limit n — > oo. Fig. 2 
clearly displays the following features. C(T) oscillates 
log-periodically around the mean value d — — ln2/lnri 
with frequency lu — — 27r/lnri. Notice that each curve 
completes about n periods (n = 8 in the figure but we 
have verified this behaviour for higher depth n as well). 
In Fig. 3a we show a typical example illustrating the de- 
pendence of the amplitudes of the specific heat on the 
scales (ri, rz)- We have plotted the maximum and mini- 
mum values of the low temperature regime of C(T) (de- 
noted respectively C + and C~), together with the mean 
value (C) = — ln2/lnri for the family r\ + T2 = 2/3. 
The oscillations decrease in amplitude as r\ decreases, 
and for t\ sufficiently small (t\ ~ 0.1) assymetries be- 
come significant. 

A point of view that allows for understanding quan- 
titatively the amplitudes of the oscillations consists in 
relating the thermodynamical properties to those of the 
spectrum. For instance, a constant value of the specific 
heat C = a is associated in general to the fact that the 
cummulative density of states (or spectral staircase) 

N(e) = I p(x)dx (12) 
Jo 

scales with energy as e a (equipartition principle) . In our 
case it can be verified that the spectral staircase grows 
approximately as e d (see Fig. 4a), and consequently the 
average specific heat is (C) = d. It is also apparent 
from Fig. 4b that the integrated density of states N(e) 
is a log-periodic function of the energy. (Similar results 
were obtained by Kimball and Frisch for the distribution 
of normal mode frequencies of fractal-based models jll] ; 
see also 0]). In fact, this features could have been antic- 
ipated by noting that N(e) also satisfies a simple scaling 
law: 

N( ri e) = N(e)/2, (13) 

whose general solution can be written as a power law 
times a log-periodic function jl3) . 

We will now show that an excellent description of the 
specific heat can be achieved by considering the first non- 
trivial correction to the bare power-law scaling for N(e): 

N(e) k e d [a + bcos(uj In £-(/))} . (14) 

In principle there is not a unique criterion for choosing 
the parameters a, b, (f>. We have determined a and b by 
requiring that the exact N(e) and its smooth approxi- 
mation (^) have the same average value and variance. 
The condition of maximum overlap between the exact 
staircase fluctuations N/e d and the cosine function fixes 
</>. 
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FIG. 3. Amplitudes of the oscillations of the specific heat 
in the log-periodic regime (T <C 1). Maximum and minimum 
values of C(T) (respectively C + and C~) for the family of 
fractal spectra n + r2 = 2/3. Also shown is the average 
specific heat (C) . For the discrete case (a) , (C) = — In 2/ In n . 
In the continuous case (b), (C) = 1 — ln(ri + r2)/lnn. 

An approximate partition function is now written in 
terms of the smoothed cummulative level density, 

/•oo 

Z(T)wj9/ e d [a + 6cos(o;lne-^)]exp(-^e)de (15) 
Jo 

= T d [aa + bb' cos(w In T — <p) — be sin(w In T - 0)] . 
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Here the constants a' ,b', d are calculated as the integrals: 

1 



a 

b' > = / dxx d e~ x I cos(cjlna;) 
' I Jo I sin(wlna;) 



(16) 



For not very small values of ri one has a' 3> b',d (e.g. 
if r x = 0.34, then a' w 0.9, 6' w 0.002, and c' w 0.0003.) 
After some straightforward manipulations, to first order 
in the small parameters b' ja! and d/a', we obtain the 
specific heat 



C(T) » d + a" cos(w In T) + b" sin(w In T) 



(17) 



This expression can be seen as a log-Fourier expansion 
of the specific heat up to second order terms. Instead 
of presenting (complicated) expressions for the constants 
a" and b" as functions of r\ and we prefer to show 



the specific heat (JIT]) for a set of selected values of r\, 
ri (Fig. 2, dotted lines). The agreement of our approx- 
imation (|17|) with the exact (numerical) calculations is 
excellent for the three upper curves and reasonably good 
for that corresponding to the smallest n (higher order 
terms might be necessary in this case). 



III. CONTINUOUS AND MULTI-SCALE 
EXTENSIONS 

For the sake of completeness let us also discuss the 
case of a spectrum constructed by iterative use of the 
rule {r\,T2) but starting from a continuous pattern as 
shown in Fig. lb. For instance, if the zero-th hierarchy 
is chosen as a continuous spectrum with uniform density 
in the interval [0, 1], then n — 1 corresponds to a spec- 
trum whose first and second bands are the intervals [0, ri] 
and [1 — 7"2, 1] respectively; and so on for increasing n. We 
have chosen the density to be uniform inside each band. 
In other words the number of states in each band is pro- 
portional to its length, whereas in the discrete case each 
"band" contains the same number of states; this is the 
essential difference between what we call continuous and 
discrete spectra. Now the partition function is written as 



An) 



(n + r 2 ) n Jo 



p{e) exp(-/3e)de (18) 



where, as in (^|), a normalization prefactor has been in- 
cluded. The analogous to (Q) is a recurrence equation for 
the density of states: 



Jn+l) 







p( n J([e-l]/ra + l) if l-r 2 <e<l 



if < e < n 
if rx < e < 1 — r2 



(19) 
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FIG. 4. (a) Integrated density of states TV (normalized to 
unity) vs. energy. The full line corresponds to the discrete 
Cantor spectrum of Fig. la with n = 0.34, r2 = 0.20, and 
n—12. The dashed line is given by e d , with d = — ln2/lnn, 
the spectral dimension, (b) Spectral fluctuations. N(e) after 
dividing by e d (full line) and the smooth approximation (dot- 



ted line) N/e = a + 6cos(o;lne - 



where u) 



-lit I In n. 



leading to the following result for the partition function 



a and b are determined by requiring that the exact and the 
smoothed fluctuations have the same average value and vari- 
ance. The condition of maximum overlap fixes <f>. 

In the low temperature regime (T <C 1), the expression 
above tends to the the scaling relation 
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ri 



n + r 2 



(21) 



Note that, as in the discrete case, the scale factor r\ is 
responsible for the period of the log-oscillations. Thus 
the discrete and continuous characteristic frequencies are 
equal. The essencial difference is the presence of ri in 
( pl| ) , which can be traced back to a different distribution 
of the spectral density (of course, when r\ — r 2 both 
cases coincide). In consequence, r 2 will also affect the 
mean value of specific heat, which is easily shown to be 
now 



(Ccont) = 1 - ha(n + r 2 )/lnn = d' 



(22) 



Equality (|22j) defines a new dimension d', which, together 
with d and df, constitute the basic set of characteris- 
tic dimensions of our problem. We remark that these 
dimensions assume different values, except for the par- 
ticular case n =■ r%. Even though the mean value ( p2] ) 
differs from its discrete counterpart (|l0|), the continu- 
ous and discrete specific heats oscillate in a similar way 
about their respective averages. However, the small-ri 
assymetries are more pronounced in the continuous case. 
These facts can be appreciated by comparing Fig. 3a and 
Fig. 3b. 

The previous analysis for the two-scale spectrum (ei- 
ther discrete or continuous) can also be generalized to the 
multiscale case. The construction of a multiscale fractal 
spectrum starts from an arbitrary discrete or continuous 
set of levels in the interval [0, 1] (n — 0). Then one makes 
M rescaled copies of the pattern n = 0, with different 
scale factors r±, . . . , r/vf. Each one of this copies is placed 
in the unit interval at positions a±, . . . , clm, so that the 
copies do not overlap (this requires ai+ri < aj+i). Itera- 
tion of this rule eventually leads to a fractal of dimension 
Si=i r i S = 1- Analogous considerations to those made 
for the two-scale case result in the following relationships 
for the discrete and the continuous multiscale cases, re- 
spectively 



1 M 

Z disc (P) = —J2e-^Z disc (0 ri ) 



(23) 



z C ant(P) = —i — J2 ne-^z^tiPn) . (24) 

These partition functions lead to the average specific 
heats (T < 1) 



(C disc ) = -InAf/lnn = d, 
(C cont ) = l-ln (X>^ / ln 



n = d' 



(25) 
(26) 



(Naturally, the d and d' we introduced in Section II are 
the M = 2 particular case of those defined above.) Once 



more the scaling exponents only depend on r\ in the dis- 
crete case and on the whole set of scaling factors rj in the 
continuous (banded) case. As in the two-scale case, these 
scale factors will be the essential ingredients for a very 
good approximate description of the thermodynamics of 
the system. 

We point out that the fractals considered in this paper 
might also be analized in their outbound and complete 
versions (in the nomenclature of ||). These variations, 
which can also be treated within our formalism, will give 
rise to a thermodynamics analogous to that described 
above. 



IV. CONCLUSIONS 

The models we have studied suggest that the hierar- 
chical organization of the energy spectra reflects itself in 
the specific heat in two ways. Simple scaling arguments 
show that the average behaviour is associated to a non- 
integer spectral dimension (d and d! in our examples), 
which in general is different from the fractal dimension 
(df). The corrections to this result are log-periodic os- 
cillations which can be traced back to the log-periodicity 
of the spectral staircase. The number of oscillations that 
can be observed is related to the hierarchical depth of the 
fractal spectrum, implying that these anomalies may ap- 
pear in systems displaying a self-similar spectrum up to 
a finite hierarchichal depth. These observations, related 
to multi-scale fractal spectra, might also be relevant in 
the case of the more realistic multifractal ones, because 
usually a few scales suffice for a good description of a 
multifractal spectrum. 

Moreover, since the effect is a consequence of the scale 
invariance of the spectrum it is quite plausible that sim- 
ilar fenomena would generically exist for bosonic and 
fermionic systems ]T^| . To support this conjecture, let 
us mention that Petri and Ruocco have observed 
fractional scaling laws when studying the (Debye) vi- 
brational specific heat of a one-dimensional hierarchical 
model. However, those authors were mainly concerned 
with mean values and did not discuss the small ampli- 
tude oscillations that can be clearly observed in their 
results. 

Even though it is not surprising that log-periodic cor- 
rections (or "complex exponents") are a natural conse- 
quence of discrete scale invariance JL5| , a contribution of 
the present paper is to have reported and analyzed exam- 
ples in which the connection between scale invariance (of 
an energy spectrum) and log-periodicity (of the specific 
heat as a function of temperature) shows up transpar- 
ently. 

Before concluding, let us comment on a possible con- 
nection of the present calculation with the recently intro- 
duced nonextensive thermostatistics |lq ]. Alemany U^A 
has suggested that this formalism could be connected to 



5 



systems with fractally structured Boltzmann-Gibbs prob- 
ability distributions. Although, for our present calcu- 
lation, we have not succeeded in making a transparent 
connection along Alemany's lines, it is worthy mention- 
ing one intriguing feature. The generalized specific heat 
C q (T) of the quantum one-dimensional harmonic oscil- 
lator does present oscillations if the entropic index 
q satisfies q < 1. In fact, C q (T)/T 1 ~ q is an oscillatory 
function of T, in a similar way C(T) is a periodic function 
of InT. 
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